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Abstract. We consider the possibility that cosmic mag- 
netic field, instead of being uniformly distributed, is 
strongly correlated with the large scale structure of the 
universe. Then, the observed rotational measure of extra- 
galactic radio sources would be caused mostly by the 
clumpy magnetic held in cosmological hiaments/sheets 
rather than by a uniform magnetic held, which was of- 
ten assumed in previous studies. As a model for the inho- 
mogeneity of the cosmological magnetic field, we adopt a 
cosmological hydrodynamic simulation, where the field is 
passively included, and can approximately represent the 
real field distribution with an arbitrary normalization for 
the field strength. Then, we derive an upper limit of the 
magnetic field strength by comparing the observed limit of 
rotational measure with the rotational measure expected 
from the magnetic field geometry in the simulated model 
universe. The resulting upper limit to the magnetic field in 
filaments and sheets is Bf s ^ 1/iG which is ~ 10 3 times 
higher than the previously quoted values. This value is 
close to, but larger than, the equipartition magnetic field 
strength in filaments and sheets. The amplification mech- 
anism of the magnetic field to the above strength is un- 
certain. The implications of such a strength of the cosmic 
magnetic field are discussed. 
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1. Introduction 

The strength and morphology of the intergalactic mag- 
netic fields remain largely unknown, because it is intrin- 



sically difficult to observe them (for recent reviews, see 
Kronberg 1994; Beck et al. 1996; Zweibel & Heiles 1997). 
While the magnetic fields inside typical galaxies are ob- 
served to be of order of 3 — 10/iG, recent observations 
indicate that the magnetic fields of the similar strength 
are also common in core regions of rich clusters (for a re- 
cent discussion see Enfilin et al. 1997). On the other hand, 
the upper limit for a large-scale field strength placed by 
the observed rotation measure (RM) of quasars is about 
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is the reversal scale of the magnetic fields in units of Mpc 
in comoving coordinates. It was assumed the fields are uni- 
form in direction within the reversal scale and varied as 
a comoving passive field in the expansion of the universe. 
The above value is close to the limit placed by the observed 
anisotropy in the cosmic microwave background radiation 
is B IGM < 6.8 x lQ- 9 (n o ft 2 ) 1 /2G (Barrow, Ferreira & Silk 
1997). Here, a uniform magnetic field was assumed, h is 
H in unit of 100 km/s/Mpc. 

According to a recent popular view based on both ob- 
servational and theoretical cosmology, the dominant non- 
linear structure in the universe is a web-like network of fil- 
aments (e.g., Bond, Kofman & Pogosyan 1996; Shectman 
et al. 1996). At a lower density contrast, however, bub- 
bly walls whose intersections are in fact filaments are the 
dominant structures. Such large scale structures form via 
gravitational instability, and then turbulent motions in- 
side the structures as well as streaming motions along the 
structures necessarily exist (e.g., Kang et al. 1994). Hence, 
magnetic fields in the early universe, if existed, should 
have been modified and amplified by those flow motions 
(e.g., Kulsrud & Anderson 1992; Kulsrud et al. 1997). As 
a result, we expect that the "geometry" of cosmic mag- 
netic field, which means the spatial distribution of the field 
strength and its orientation, should be correlated with the 
large-scale nonlinear structures of the universe. In other 



2 



Cosmic Magnetic Fields 



words, the field strength increases with the matter den- 
sity and its orientation tends to align with the sheets and 
filaments, while its random component is associated with 
the local turbulent motions. The magnetic fields should 
be strong along the walls of cosmic bubbles, even stronger 
along the filamentary superclusters, and the strongest in- 
side the clusters of the galaxies, while they are very weak 
inside the voids. 

Based on the proposition that the large scale magnetic 
field is correlated with the large scale structure of the uni- 
verse, we consider a way to estimate the field strength 
which relies on the observational data. 

The cosmic magnetic field together with free elec- 
trons in the intergalactic medium (IGM) induces the Fara- 
day rotation in polarized radio waves from extra-galactic 
sources. The observational RM data of quasars show a 
systematic growth of RM with rcdshift, z, and limit RM 
to ~ 5 rad mT 2 or less at z = 2.5 (Kronberg & Simard- 
Normandin f976; Kronberg f994 and references therein). 
Adopting a model for the distribution of the large scale 
magnetic field, these data can be used to constrain the 
strength of the magnetic field. For example, the upper 
limit estimated by Kronberg (1994) was based on a model 
in which the field is uniform within the reversal scale, the 
orientation is random, and there are N = l s /L rev rever- 
sals along the line of sight to a source l s Mpc away from 
us. 

In the present paper, we re-derive the upper limit by 
taking a new model in which the intergalactic magnetic 
field is mostly confined within filaments and sheets but 
very weak inside voids, rather than the simple model of 
uniform field. Of course the field would be strongest in- 
side clusters, but their contribution is usually excluded in 
the observed RM. If we take only the geometric consid- 
eration that the field distribution along the lines of sight 
to radio sources has a small filling factor, high only in- 
side filaments and sheets but low inside voids, then it is 
obvious that the expected magnetic field strength in fila- 
ments and sheets should be higher than the value derived 
from the uniform field model. But this geometric model is 
too simple, since the electron distribution as well as the 
field direction are uncertain. Here, we take a more prac- 
tical approach by adopting numerical data in a simulated 
universe which includes a magnetic field, which is evolved 
during the large scale structure formation. With the mag- 
netic field distribution, specific to this simulation model, 
we argue that the magnetic field strength in filaments and 
sheets is limited to Bj s ^ 1 ^G. We should emphasize 
that this limit should be approximately valid regardless of 
details of any model for the origin of the large scale mag- 
netic field, provided that it is correlated with the large scale 
structure of the universe. 

In the next section, we describe the model simulation 
and the resulting magnetic field geometry. In §3, we ex- 
plain the procedure to calculate the upper limit of the 
magnetic field strength in filaments and sheets constrained 



by the observed RM. In the final section, we discuss the 
implications of the possible existence of ~ 1 ^G or less 
magnetic field in filaments and sheets. 

2. Model Simulation 

In this section we describe one particular simulation of 
the origin and growth of magnetic fields, but will retain 
at the end only the geometry and the relative strength of 
the magnetic fields across the cosmological structures. It is 
this structure for which we require a realistic simulation 
to probe the universe for its Rotation Measure across a 
large number of simulated paths. 

The generation of the large scale magnetic field by 
the Bicrmann battery mechanism (L. Biermann 1950) and 
its subsequent evolution has been followed by solving the 
equation 

<9B „ , „. cVp e x Vn e , . 

along with the cosmological hydrodynamic equations. 
Here, p e is the electron pressure and n e is the electron 
number density. The simulation has been done with a ver- 
sion of the cosmological hydrodynamic code which extends 
the one described in Ryu et al. (1993) by incorporating the 
above equation "passively". That is, the Lorenz force term 
in the momentum equation has been ignored. Alterna- 
tively, the equation for the magnetic field could have been 
incorporated by solving the full MHD equations. Then the 
code would have been more diffusive. 

A standard cold dark matter (CDM) model universe 
with total fl = 1 has been simulated in a periodic box with 
(32/i _1 Mpc) 3 volume using 128 3 cells and 64 3 particles 
from Zi = 20 to Zf = 0. The values of other parameters 
used are fib = 0.06, h = 1/2, and erg = 1.05. A more 
detailed description of the simulation can be found in in 
Kulsrud et al. (1997). 

In the simulations which include the evolution of bary- 
onic matter as well as that of dark matter, accretion 
shocks form in the infalling flows towards the growing 
nonlinear structures such as sheets, filaments and clus- 
ters (e.g., Kang et al. 1994; Ryu & Kang 1997b). We note 
that some evidence has been found in radio relic sources, 
that these accretion shocks have actually been detected 
(Enfilin et al. 1998). The properties of the shocks depend 
upon the power spectrum of the initial perturbations on 
a given scale as well as the background expansion in a 
given cosmological model. These shocks are the sites of 
the generation of weak seed field by the battery mecha- 
nism represented by the last term in the right hand side 
of Eq. (1). The evolution of the strength of the seed field 
generated was plotted in Figure 1 of Kulsrud et al. (1997). 
It grows monotonically to the mass averaged value of or- 
der 10 _21 G by z <~ 2 — 3, and then levels off without 
further increase. The saturation is believed to be due to 
the finite numerical resistivity inherent in the numerical 
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scheme used to solve Eq. (1). At the saturation, the mag- 
netic energy is much smaller (by ~ 10 30 ) than the kinetic 
energy. So the dynamical influence of magnetic field into 
flow motion can be ignored, and neglecting the Lorenz 
force in the momentum equation is justified. 

To illustrate the general characteristics of the field ge- 
ometry in relation to the matter distribution and flow ve- 
locity, a two-dimensional cut of the simulated box is pre- 
sented in Figure f. The plot shows a region of 32/i _1 x 
20/i- 1 (Mpc) 2 with a thickness of 0.25h _1 Mpc at z = 0. 
The three panels display baryonic density contours, veloc- 
ity vectors, and magnetic field vectors. 

We can see accretion shocks around the high density 
structures of clusters, filaments, and sheets in the density 
contours. The locations of the shocks can be seen more 
clearly in the velocity vector plot where accretion motions 
stop suddenly In the highest density regions of clusters, 
flow motion is turbulent. Also we can see the streaming 
flow motion along the structures, that is, the accretion mo- 
tion towards clusters in filaments and the accretion motion 
towards filaments in sheets. On the top of the streaming 
motion, we expect some amount of turbulent motion which 
is induced by the cooling, gravitational, and nonlinear thin 
shell instabilities in filaments and sheets (Vishniac 1994; 
Anninos, Norman & Anninos 1995; Valinia et al. 1997). So 
the magnetic field in filaments and sheets could have been 
stretched by the streaming motion as well as amplified by 
the turbulent motion. As a result, the amplified magnetic 
field in filaments and sheets is expected to be aligned with 
the structures, as shown in Figure 1. 

Here, we propose that the geometry of the magnetic 
field in the simulated universe, with an arbitrary normal- 
ization for the field strength, represents the geometry of the 
real magnetic field correlated with the large scale structure 
of the universe. Thus we take from the simulation only the 
geometry of the field distribution and its relative strength 
across the spatial distribution. 



3. Upper Limit on the Magnetic Field in Filaments 
and Sheets 

We have calculated the expected RMs induced by the mag- 
netic field in filaments and sheets in the simulated universe 
and compared them with the observed upper limit of RM 
of distant quasars. 

Due to the field reversal of randomly oriented magnetic 
field along lines of sight, the RMs from sources at a given 
rcdshift have a Gaussian distribution centered at zero and 
its standard deviation gives a statistical measure for the 
field strength. The standard deviation of the distribution 
of the RMs from quasars at z ~ 2.5 is RM <, 5 rad mT 2 
(Kronberg 1994). 
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The RM which measures the amount of the rotation is 
given by 

where ip is the rotated angle, X(l) is the wavelength of 
the polarized wave along the propagating path, I, A is the 
observed wavelength, and B\\ is the magnetic field com- 
ponent along I (e.g., Shu 1991). The RM for a source at 
a given redshift z can be calculated numerically by the 
following line integral 

RM (jj^j = 8Axl ° 5 J a ' h eB\\dl = 9 - 2 xh 2 J^ ' Sl b B\\dl (3) 

where I is the comoving length in the unit of Mpc, B\\ is the 
proper magnetic field in [iG, n e is the comoving electron 
number density in cm~ 3 , is the baryonic density as 
a fraction of the critical density, and x 1S the ionization 
fraction. The integration length I is related to the redshift 
by the following equation 

rf/ = --^(l + z)- 1 (l + r! z)- 1 / 2 dz. (4) 

In the calculation, we have assumed the gas inside fila- 
ments and sheets is fully ionized (i.e., \ = !)■ The simu- 
lated distributions of the gas and magnetic field at z = 0, 
0.5, 1, and 2 have been used, and they have been inter- 
polated in between to integrate the above equation up to 
z = 2.5. 

For the relative growth of the field strength for 2.5 > 
z > 0, however, we need to adopt a model, since the nu- 
merical simulation could not follow the growth of mag- 
netic field in detail due to numerical resistivity. Again, 
the normalization of the field strength is arbitrary. Con- 
sidering that theoretical studies on dynamo processes on 
large scales do not yet provide a unambiguous model, we 
assume for simplicity that the magnetic field has reached 
an "equipartition" with the energy in turbulence. Then its 
energy should be proportional to the kinetic energy which 
is, in turn, proportional to the gravitational potential en- 
ergy. Then, the self-similar solution for an £1 = 1 universe 
gives 

- 2(2rce-l) 

Bcx(l + z)^^, (5) 

where n = 1 is for sheets, 2 for filaments, and 3 for clus- 
ters (Fillmore & Goldreich 1984). Here, e is the parameter 
which characterizes the initial perturbation of the self- 
similar solution and has < e < 1. Hence, for filaments 
and sheets, we expect 

% s oc(l + z)« (6) 

with ^ q 1.5. However, we have also considered the 
model in which the magnetic field grows as time passes 
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Fig. 1. Two-dimensional cut of the simulated universe at z = 0. The plot shows a region of 32/i _1 x 20/i _1 Mpc 2 with a thickness 
of 0.25/i _1 Mpc, although the simulation was done in a box of (32/i -1 Mpc) 3 volume. The first panel shows baryonic density 
contours, the second panel shows velocity vectors, and the third panel shows magnetic field vectors. In the first panel, the solid 
lines contour the regions of pb > pb and the dotted lines contour those of 10 _1,2 p6 < pb < pb- In the third panel, the vector 
length is proportional to the log of magnetic field strength. 
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(Welter, Perry, & Kronberg 1984). So, we have explored 
the cases of —1.5 < q < 1.5. This large range of the param- 
eter q explored here should encompass many models for 
the origin of magnetic fields in the cosmos (See, e.g., the 
discussions in L. Biermann 1950; Parker 1958; Rees 1987; 
Beck et al. 1996; P. Biermann 1997; Kronberg & Lesch 
1997). 

With the above prescription for the evolution of mag- 
netic field in filaments and sheets, we have integrated Eq. 
(3) up to z = 2.5 along 3 x 10 4 randomly chosen paths. 
Figure 2 shows the distribution of the expected RMs for 
the model with q = 0. The magnetic field is normalized so 
that the standard deviation of the Gaussian fit (the solid 
line) to the expected RMs matches the observational up- 
per limit of RM = 5 rad m -2 . With this normalization for 
the field strength, the upper limits in the rms of magnetic 
field in filaments and sheets can be calculated for several 
models. They are given in Table 1 for the several values 
of q, and they are ~ 1 fiG for h = 1/2. 




-20 -10 10 20 
RM(rad/m 2 ) 

Fig. 2. Histogram of the fraction of the expected RM along 
3 x 10 4 randomly chosen paths up to z = 2.5 in the simulated 
universe. The solid line represents the Gaussian fit. The ex- 
ponent in the temporal evolution of the rms of magnetic field 
in filaments and sheets, q, is assumed to be zero. The mag- 
netic field strength is normalized so that the standard devia- 
tion of the Gaussian fit matches the observational upper limit 
of RM = 5 rad m~ 2 . 



Table 1. Upper Limits on the Magnetic Field Strength in Fil- 
aments and Sheets 

q a Bz =0 .j s (jj,G)hl 5 b 

1.5 0.27 

1 0.42 

0.5 0.63 

0.89 

-0.5 1.15 

-1 1.39 

-1.5 1.58 

a The exponent in the temporal evolution of magnetic field in 
filaments and sheets, B z= q.j s oc (1 + z) q . 

b The rms of magnetic field in filaments and sheets at the 
present epoch (z = 0) in units of )j,G for h = 1/2. 



10 

1 

^1 0.1 
a 0.01 
0.001 

: 



E 1 1 


1 




i i 




iii = 

z=o : s 

Mil 


1 1 


1 


1 1 
, , , 


1 1 
i i 


\ 

1 


1 1 1 



8 16 24 32 

l(h- l Mpc) 

Fig. 3. Distributions of baryonic density and magnetic field 
strength at z = along the line drawn in the first panel of 
Figure 1, when the exponent in its temporal evolution is, q, 
zero. The density is given in units of the critical density, while 
the magnetic field is given in units of [iG for h — 1/2. 

These numbers clearly demonstrate that the essential 
result is nearly independent of any of the details about 
the origin of the magnetic field, since for the entire range 
of values tested in q we obtain to within a factor of a few 
the same result: The upper limit to the magnetic field in 
the sheet is of order Bf s 55 1/xG. 



Figure 3 shows the distributions of gas density and 
magnetic field strength along the line drawn in the first 
panel of Figure 1 for the case with q = 0. We can see 
the field strength correlates well with the gas density and 
ranges 0.2 < B < 2fiGh^ 2 5 in the regions of filaments and 
sheets. 



4. Conclusion and Discussion 

Our result indicates that, with the present value of the ob- 
served limit in RM, the existence of magnetic field of up to 
~ 1 fiG in filaments and sheets can not be ruled out, if the 
cosmic magnetic field is preferentially distributed in these 
nonlinear structures, rather than uniformly distributed in 
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the intergalactic space. We emphasize that this result is 
independent of the details of how the magnetic field origi- 
nated. However, it is not certain if the dynamo processes 
such as the one considered in Kulsrud et al. (1997) can 
amplify a week seed field to the above field strength. 

Interestingly, there is confirmation of such a magnetic 
field in one case. Kim et al. (1989) reported the possible 
existence of a intercluster magnetic field of 0.3 — 0.6 /iG 
in the plane of the Coma/Abell 1367 supercluster. They 
reached the result from the observation using high dy- 
namic range synchrotron images at 327 MHz. Clearly, this 
result should be confirmed in other regions of the sky with 
the newly available possibilities of low radio frequency 
intcrferometry such as with the Giant Meterwave Radio 
Telescope (GMRT) in India, or the new receiver systems 
on the Very Large Array (VLA) in the USA. 

We may compare the above upper limit with the 
strength of the magnetic field whose energy is in equipar- 
tition with the thermal energy of the gas in filaments and 
sheets. The equipartition magnetic field strength can be 
written as 

The fiducial values T = 3 x 10 6 K and = 0.3p c may 
considered to be appropriate for filaments (see, Kang et 
al. 1994 and Figure 3). The values appropriate for sheets 
should be somewhat smaller. So, our upper limit set by 
RM is close to, but several times larger than, the equipar- 
tition magnetic field strength. 

Filaments and sheets with the above fiducial temper- 
ature and gas density can contribute radiation at X-ray 
wavelengths to the cosmological background. Using the 
above values and a typical size of 10/i _1 Mpc, the lu- 
minosity of filaments and sheets in the soft X-ray band 
(0.5— 2kcV) is expected to be ^ 3 x 10 41 erg s _1 . Actually, 
Soltan et al. (1996) and Miyaji et al. (1996) reported the 
detection of extended X-ray emission from structures of a 
comparable size and a luminosity of ~ 2.5 x 10 43 erg s _1 , 
which is correlated to Abell clusters. These observations 
suggest that the temperature and gas density outside clus- 
ters could even be considerably larger than the fiducial 
values assumed above. In any case, if these observations 
can be further confirmed, and the equipartition of mag- 
netic field and gas energies is assumed, it would imply the 
existence of a <~ lp,G magnetic field on large scales outside 
galaxy clusters. 

The possible existence of strong magnetic field of 1 /iG 
or less in filaments and sheets has many astrophysical im- 
plications, some of which we briefly outline in the follow- 
ing: 

The large-scale accretion shocks where seed mag- 
netic fields could be generated are probably the biggest 
shocks in the universe with a typical size ^ a few (1- 
10) Mpc and very strong with a typical accretion ve- 
locity ^ a few 1000 km s _1 . The accretion velocity onto 



the shocks around clusters of a given temperature, or a 
give mass to radius ratio M c i/R c i, is smaller in a uni- 
verse with smaller o , and is given as v acc « 0.9 — 1.1 x 
10 3 kms- 1 [(M c/ /i? c/ )/(4 x lO^Mo/Mpc)] 1 / 2 in model 
universes with 0.1 < fl < 1 (Ryu & Kang 1997a). With 
up to ~ 1 pG or less magnetic field around them, the 
large-scale accretion shocks could serve as possible sites 
for the acceleration of high energy cosmic rays by the 
first-order Fermi process (Kang, Ryu & Jones 1996; Kang, 
Rachen & Biermann 1997). Although the shocks around 
clusters would be the most efficient sites for acceleration, 
those around filaments and sheets could make a signif- 
icant contribution as well (Norman, Melrose & Achter- 
berg 1995). With the particle diffusion model in quasi- 
perpendicular shocks (Jokipii 1987), the observed cosmic 
ray spectrum near 10 19 eV could be explained with rea- 
sonable parameters if about 10 -4 of the infalling kinetic 
energy can be injected into the intergalactic space as the 
high energy particles (if an E~ 2 spectrum of cosmic rays 
is assumed; for a slightly steeper spectrum as suggested by 
radio relic sources, the efficiency is closer to 0.1, EnBlin et 
al. 1998). 

The discoveries of several reliable events of high en- 
ergy cosmic rays at an energy above 10 20 eV raise ques- 
tions about their origin and path in the universe (a re- 
cent review is P. Biermann 1997), since their interaction 
with the cosmic microwave background radiation limits 
the distances to their sources to less than 100 Mpc, per- 
haps within our Local Supercluster. The Haverah Park 
and Akeno data indicate that their arrival directions are 
in some degree correlated with the direction of the Super- 
galactic plane (Stanev et al. 1995; Hayashida et al. 1996; 
Uchihori et al. 1996). In Biermann, Kang & Ryu (1996), 
we noted that if the magnetic field of ~ 1/iG or less ex- 
ists inside our Local Supercluster and there exist accre- 
tion flows infalling toward the supergalactic plane, it is 
possible that the high energy cosmic rays above the so- 
called GZK cutoff (£>5x 10 19 eV) can be confined to 
the supergalactic plane sheet, an effect analogous to so- 
lar wind modulation. In each case a shock front pushes 
energetic particles upstream as seen from its flow. Obvi- 
ously, this effect would occur only for a small part of phase 
space, namely those particles with a sufficiently small ini- 
tial momentum transverse to the sheet. This would explain 
naturally the correlation between the arrival direction of 
the high energy cosmic rays and the supergalactic plane. 
Also, confinement means that for all the particles cap- 
tured into the sheets, the dilution with distance d is 1/d 
instead of 1/d , increasing the cosmic ray flux from any 
source appreciably with respect to the three-dimensional 
dilution. So we may see sources to much larger distances 
than expected so far. On the other hand, particles with 
a larger initial momentum transverse to the sheet would 
be strongly scattered, obliterating all source information 
from their arrival direction at Earth. 
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With the magnetic field in the intergalactic medium, 
charged particles passing it would not only experience de- 
flection. It would also smear out their arrival direction as 
well as delay their arrival time. Plaga (1995) suggested 
that an exhibition of this latter effect would be the delay 
of the arrival times of 7-rays from a cascade caused by 
photons from highly time-variable extragalactic sources. 
However, with strong fields of ^ 1 /iG or less in filaments 
and sheets intervening between the sources and us, the 
expected consequence would be a strong smearing of the 
sources rather than the delay of arrival times (Kronberg 
1995). But, for details, calculations following the propa- 
gation of photons and charged particles should be done. 

Another interesting exploration is to study the radi- 
ation emission arising from filaments and sheets strong 
magnetic field of 1 fiG or less. We noted that the orig- 
inal estimate of the magnetic field in the plane of the 
Coma/Abell 1367 supercluster was based on a synchrotron 
radio continuum measurement (Kim et al. 1989). We also 
noted that the thermal Bremsstrahlung emission in the 
soft X-ray band may provide an interesting check on the 
work presented here. Whether these structures can con- 
tribute at other wavelengths, such as 7-ray energies, to 
the cosmological background remains an an unanswered 
but challenging question at this time. 

These issues will be discussed elsewhere. 
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